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We consider a long superconductor-ferromagnet-superconductor junction with one spin-active 
region. It is shown that an odd number of Cooper pairs cannot have a long-range propagation when 
there is only one spin-active region. When temperature is much lower than the Thouless energy, the 
coherent transport of two Cooper pairs becomes dominant process and the superharmonic current- 
phase relation is obtained (7 oc sin 20). 



The interplay between superconducting and ferromag- 
netic ordering has been the subject of intensive theoreti- 
cal and experimental research [l|-|9| . It has been predicted 
that the hybrid systems containing superconductors (S) 
and ferromagnets (F) allow the realization of the Joseph- 
son 7r-junctions, / oc sin(</> + tt) The spectral decom- 
position of Josephson current-phase relation (CPR) gives 
I = h sin <p + 1-2 sin 2<fi ... . At the transition between the 
and the n phases, I\ vanishes, and it is possible to 
obtain dominant second harmonic (I2) in the CPR. Un- 
fortunately, the 0-7T coexistence is very sensitive to the 
temperature changes and interface roughness ■ 

The experimental realizations of 7r-junctions remained 
elusive for a long time; the breakthrough came with the 
fabrication of weak ferromagnets Indeed, it has been 
recognized that the proximity effect in a fcrromagnct 
is short-ranged. The electron and hole excitations ac- 
quire a nonzero relative phase in the ferromagnet between 
the scatterings from the two superconductor-ferromagnet 
(SF) interfaces 0; the different orbital modes acquire dif- 
ferent phases that add up destructively after summation. 
A detailed analysis demonstrates that the proximity 
effect in a superconductor-ferromagnet-superconductor 
(SFS) junction is suppressed algebraically in the ballistic 
regime, and exponentially in the diffusive one 

Quite recently, it was proposed that a SFS junction 
with an inhomogeneous magnetization in the F layer can 
generate a triplet pairing and can support a long-range 
Josephson current [J] . Subsequent theoretical and exper- 
imental research showed that in order to have dominant 
triplet pairing a SFS junction with two spin-active re- 
gions is required 

Motivation for this work is our previous numerical cal- 
culation , where clean and moderately disordered SFS 
junctions were considered (with a one spin-active inter- 
face), and the dominant second harmonic was obtained. 

In this work, we consider a long SFS junction in the 
ballistic regime. Assuming the presence of a spin-active 
region on only one SF interface, we show that only the 
phase coherent transport of an even number of Cooper 
pairs is not suppressed by the exchange field. In par- 
ticular, the dominant contribution to the Josephson cur- 
rent stems from the transport of two Cooper pairs. As 
a consequence, the two times smaller flux quantum is 



obtained, leading to more sensitive quantum interferom- 
eters (SQUIDs) [ll[ and the half-integer Shapiro steps 
that can be experimentally observed [9j. Another inter- 
esting property is the coexistence of integer and half- 
integer fluxoid configurations in SQUIDs, corresponding 
to the minima of the triple- well potential energy [Tlj : 
this can be potentially useful for experimental study of 
the quantum superposition of macroscopically distinct 
states [u|. Also, junctions with a nonsinusoidal current- 
phase relation are shown to be promising for realization of 
"silent" phase qubits Last but not least, this result 
enables robust realization of so-called 95-junctions [f"| . 

It should be stressed that in contrast to the case of 
the 0-7T transition the discovered effect is very robust: 
it is insensitive to a weak disorder [Ioj [. temperature 
changes, and the interface roughness. Nevertheless, rela- 
tively transparent interfaces arc required in order to ob- 
serve the effect. 

Before we proceed with a quantitative analysis of the 
aforesaid effect, let us first give a simple and intuitive de- 
scription. Note that when a SF interface is spin-active, 
there are two possibilities for Andreev reflection: the nor- 
mal Andreev reflection (the spin projections of an elec- 
tron and the reflected hole are opposite) and the anoma- 
lous Andreev reflection (an electron and the reflected hole 
have the same spin projections) [l5|. We consider sepa- 
rately the phase coherent transport of one (Ji) and two 
(J2) Cooper pairs. The transport of a single Cooper pair 
is suppressed by the exchange field because the electron 
and the Andreev reflected hole have opposite spin pro- 
jections (see Fig. [I]). On the other hand, the transport of 
two Cooper pairs has a long-range contribution stemming 
from two normal and two anomalous Andreev reflections. 

We consider a simple model of a ballistic SFS junc- 
tion consisting of two conventional (s-wave) supercon- 
ductors, a uniform single-domain ferromagnet and only 
one spin-active region. Andreev reflection requires rel- 
atively transparent SF interfaces, thus for simplicity we 
assume them to be fully transparent. The spin-active 
region consists of a ferromagnetic spacer layer with the 
magnetization noncollinear to that of the F layer. 

The Josephson current is calculated using the scatter- 
ing approach [l6[ . The knowledge of scattering matrices 
(S-matrices) of both SF interfaces is sufficient to obtain 
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FIG. 1. (color online) The first two harmonics in the Joseph- 
son current-phase relation. The first one (7i) consists of two 
normal Andreev reflections from SF interfaces. The second 
one (72) has two contributions: the one with four normal An- 
dreev reflections (short-range), the other with two normal and 
two anomalous Andreev reflections (long-range, total phase 
acquired in the F layer is zero). The solid (dotted) lines rep- 
resent electron (hole) excitations in F layer. The red (vertical) 
arrows represent spin projections, while the black (horizontal) 
arrows denote the excitation velocity direction. The right SF 
boundary (hatched) is spin-active. 



the Josephson current in the ballistic regime. Each of 
these matrices relates the amplitudes of the excitations 
propagating towards the corresponding SF interface to 
the excitation propagation away from it. There is no or- 
bital channel mixing in the ballistic regime, but only mix- 
ing of different spin channels (due to spin-active region) 
and the particle-hole mixing (due to superconductors). 
Ergo, the dimension of the S-matrix is 4 x 4 — we write 
all matrices in the Kronecker product of particle-hole and 
spin spaces. The Josephson current is given by fl6| 
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J2 lndet[l - K{iu n )n' (iuj n )l (1) 

where TV and 1Z are S-matrices of the two SF interfaces, 
while the phases acquired upon propagation through the 
F region are included in one of these matrices. uj n = 
(2n + l)7rfcsT are the Matsubara frequencies and <f> is 
the phase difference between the two superconductors. 

In order to calculate the S-matrix, one has to solve 
the Bogoliubov-de Gennes equation for each SF inter- 
face. Here we take a simpler approach and express the S- 
matrix (TV) in terms of a S-matrix of a SN interface (3.[l5j. 
where N stands for normal-nonferromagnetic layer. By so 
doing, we neglect the difference in the number of spin-up 
or -down modes in the ferromagnet. This approximation 
is also assumed in the quasiclassical approach and isjus- 
tified for a weak exchange field in the ferromagnet [17J . 

The S-matrix of a transparent SN interface reads 



TZ A (e,<p) = a(e) 







MJ2<1 





(2) 



where a(e) = e - J arccos ( £ / A ) an( j a2 [ s the second Pauli 
matrix. For the ferromagnetic spacer layer (spin-active 
region), the S-matrix is 



Sf = 



u 
u 



(3) 



where m = (sin 9 cos i[>, sin 9 sin ip, cos 9) is the magneti- 
zation orientation in the ferromagnetic spacer layer. The 
difference of the phase shifts of spin-up and spin-down 
electrons upon propagating through the ferromagnetic 
spacer layer is denoted by p = — , while 77 = iaj- + . 
Here the spin-up (spin-down) is defined with respect to 
the magnetization axis in the F layer. These phases de- 
pend on the orbital channel index /1, but for the sake of 
notational simplicity we have suppressed the index. We 
also assume that the ferromagnetic spacer layer thick- 
ness is much smaller than the superconducting coherence 
length (V <C £s)j so that the energy dependence of p and 
r\ can be ignored. Combining these scattering matrices 
we obtain the scattering matrix of the SF interface with 
the spin-active region 
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For the SF interface without the spin-active region, we 
also include the phases acquired in the F layer and obtain 



with T 



K'(e)=T{e)K A (e,0)T( S ), 
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the longitudinal component of wavevector in the orbital 
mode p (spin- up/down), and L is the F layer thickness. 

In order to perform the integration over orbital modes, 
we put p M = Z'/cos0, where Z' = 2h'L' /(Hvp); h' is the 
exchange energy in the spacer layer which is assumed 
to be small (h! -c Ef), and is the angle between 
the excitation velocity and the junction axis. Also, in 
the Andreev approximation k^^^(e) = (hp ± Ii/Hvf + 
e/Hvf)/ cos 0, where h is the exchange energy in the F 
layer (h <C Ef)- Setting L' = (no spin-active layer) 
recovers the result previously obtained from the quasi- 



classical approach where the following expression 
relates the scattering to the quasiclassical approach (see 
also Ref. EF 
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where Rn is the normal resistance and g a is the normal 
Green function in the ferromagnet fill ]. 

In the general case — with one spin-active region — upon 
inserting Eqs. (fj|, ^ into Eq. ([T]) the formula for the 
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Joscphson current is obtained. We introduce a new vari- 
able u = arcsinh(w/A), and integrate over orbital modes 

,7r/2 

1=—^- > ' dQ sin9 cos9 (8) 
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with x = arccos(Re[r-^e- 4Z / cose ]) and Z = 2hL/(h 



v F ) 



denotes the element (2, 1) of the matrix in Eq. ©. 
We now consider a long SFS junction, and show that 
only even harmonics in the CPR arc long-range. In this 
case (L ^> £s), the first term in the argument of the hy- 
perbolic tangent (u) can be neglected. At zero tempera- 
ture, the summation over ui n can be replaced by an inte- 
gration; the expression for the Josephson current reads 



I = 



f \ -> 



eR N L ^ 



(-l) fc 4sin(fc</>). 



(9) 



In the last equation, the spectral weights Ik quantify the 
contribution to the current coming from the phase coher- 
ent transport of k Cooper pairs across the barrier 



h = - 
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where Tk is the Chebyshev polynomial of the first kind 
and x = 1/cose. Wc have used the identity Imtanh z = 
2 E£liM) fc ± me ~ 2fez for obtaining Eqs. ([MTU]) [lj|. 

In order to avoid cumbersome expressions, we concen- 
trate on the case with mutually orthogonal magnetiza- 
tions in the F and the ferromagnetic spacer layer (8 — 
7r/2). Later we will show that all our conclusions are valid 



for arbitrary (but not too small) 9. Now, 



ae 



cos(Z'x) 



and the expression for the spectral weights reads 



h = 



dx 
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EE 



(11) 



n=0 m— 

x [cos(Zx) cos(Z'x)} k - 2in - rn \ 

where [^J denote the largest integer not greater than x. 

We perform the integration in Eq. (fTTj) by expanding 
the integrand as a sum of cosines. When k is odd, every 
term in that sum depends on Z\ for a long SFS junction 
and for a reasonably strong exchange energy (h > A) 
we have Z S> 1 and J ± cos{Zx)/x A dx = — sin Z/Z + 
0(1/Z 2 ). Hence, odd harmonics are suppressed by the 
factor 1/Z. On the other hand, for even k we find terms 
that are independent of Z (the exchange energy in the F 
layer). Thus wc write I 2k = I^k + °\ l l z )- Thc lon g" 
range component of even harmonics {I^) 1S given by 
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FIG. 2. The long-range component of the second (thick 
solid curve) and the fourth (thin solid curve) harmonic of 
the Josephson current-phase relation. For Z' S> 1, both the 
curves converge to constant values (dotted lines), I2 and I4. 



This result shows that even harmonics dominate in a 
long junction. It should be noted that even harmonics 
can dominate for any junction length if the parameters 
of the spin-active region are chosen in such a way that 
normal Andrccv reflection vanishes on one SF interface 
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0). This can be directly seen from Eqs. ([ MTU]) . 



Figure [2] shows the dependence of the first two long- 
range harmonics (I2 and I4) on Z'. For large values 
of Z' both curves converge to the constant values — 
l-i and J4, respectively. Hence, we see that for a cer- 
tain range of the ferromagnetic spacer layer parameters 
[1 <C L'b! I (Hvf) *C h'/ A] the Josephson current reads 

I=^zZ^ R ^^) + 0(—, — ), (13) 
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where li? arc independent of thc junction parameters 
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Equations (|13H14p assert that for a long SFS junction 
with one spin-active region, the long-range part of the 
Josephson current depends only on thc Thoulcss energy 
Consequently, thc long-range part of the current is not 
suppressed by fluctuations of the ferromagnetic barrier 
thickness (interface roughness). The current- phase rela- 
tion is depicted in the inset of Fig. [3] 

In the general case, when the value of Z' is arbitrary 
(but L' <^ £5), thc supcrcurrent dependence on Z' is 
given by Eq. (|12p . Again, even harmonics are dominant. 
The free energy of Josephson junction is given by F(<j)) oc 
$q I((j>)d4>; the ground state of the junction is degenerate: 
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different. In the case of two spin-active interfaces, all 
harmonics in the Joscphson current-phase relation are 
long-range (and the first one is dominant) while in 
the case of one spin-active interface, we find that only 
even harmonics are long-range (and the second one is 
dominant). Some repercussions of the discovered effect 
are: half- integer Shapiro steps Q , the coexistance of in- 
teger and half- integer fluxoid SQUID configurations [11 1 
and robust realization of the (^-junctions 14 1. 

I acknowledge C. Bruder and C.W.J Beenakker for 
reading the manuscript. I am especially grateful to 
Z. Radovic and V. M. Stojanovic for many useful discus- 
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FIG. 3. (color online) The spectral weights of the current- 
phase relation for Z' = 0.5 and Z — 1000, and for three 
values of the relative angle between magnetizations: 6 = 
7r/10, it/4, 7t/2. The inset depicts the generic CPR given 
by Eq. (|13[1 . The Josephson current is normalized by 
e-Rjv/(4-E-r), where Et = Hvf/L is the Thouless energy. 



the and w states have the same energy. The degeneracy 
is lifted only by algebraically small factors (~ 1/Z). 

The dependence of the spectral weights on the misori- 
entation angle is depicted for three values of 9 in Fig. 
[3l We observe that even harmonics dominate also for 
6 < ir/2. As the relative angle between magnetizations 
approaches zero, the amplitudes of even harmonics lessen 
and eventually the odd and even harmonics become com- 
parable. Therefore, one can tune the ratio I2/I1 by 
changing the angle 9, while the sign of I\ can be changed 
by adjusting the thickness of the F layer Q. This en- 
ables robust realization of ^-junctions [Tij . contrary to 
previously discussed realizations with 0-tt transition. 

At finite temperatures a new length scale appears — the 
normal metal coherence length £jv = hvp / (2ttT). For 
high temperatures, T 3> Hvf/L = Et (i.e. L 3> £,n), 
only the first term (n = 0) in the summation over Mat- 
subara frequencies contributes to the current given by 
Eq. ([8]). After performing the integration over orbital 
modes, assuming /i>T, we obtain 



8ttT 
eR 



N 



sin(Z - Z') 
Z 



sin (Z + Z') 9 9 

— - sir - 

Z 2 



(ttT + VA 2 +7T 2 T 2 ) 2 



-e ~ L ^ N sine 



(15) 



We conclude that the first harmonic always dominate in 
the /wg/i-tcmperaturc limit, because the higher harmonics 
are suppressed by the factor e~ k ^ N ^ L (k > 1) [2(|; in this 
case the supercurrent has only the short-range part 

In conclusion, we have shown that SFS junctions with 
one and two spin-active SF interfaces are qualitatively 
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